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Hamiltonian operator for loop quantum gravity coupled to a scalar field 

Emanuele Alescij^] Mehdi Assanioussif] Jerzy Lewandowskij^] and Ilkka MakinenJ^] 
Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland 


We present the construction of a physical Hamiltonian operator in the deparametrized 
model of loop quantum gravity coupled to a free scalar field. This construction is based on 
the use of the recently introduced curvature operator, and on the idea of so-called special 
loops. We discuss in detail the regularization procedure and the assignment of the loops, 
along with the properties of the resulting operator. We compute the action of the squared 
Hamiltonian operator on spin network states, and close with some comments and outlooks. 


I. INTRODUCTION 

General relativity in Ashtekar-Barbero variables mm can be cast in an SU(2) Yang-Mills 
theory and treated as a Hamiltonian system with constraints consisting of the Gauss (gauge con¬ 
straints), spatial diffeomorphism and Hamiltonian constraints. Canonical loop quantum gravity 
GHE], which is an attempt of quantization a la Dirac [7j of general relativity, has successfully 
completed the construction of a kinematical Hilbert space and the implementation of the Gauss 
constraints and the spatial diffeomorphism constraints [8] in the quantum theory, leading to 
a gauge and spatial diffeomorphism invariant Hilbert space Tf^g- The treatment of the last 
constraints is a more complicated task. The Hamiltonian has been regularized and promoted to 
an operator acting on Tf^g by Thiemann |9j improving earlier attempts iTIj. however even if 
the general structure of the solutions to the Hamiltonian constraints is known, it is very difficult 
to define the physical Hilbert space. The issues are conceptual and technical. 

Conceptual, because the Hamiltonian is not preserving H^g and even if attempts to deal 
with the absence of a physical Hilbert space have been explored HD, this problem has led to 
new research directions, in particular the master constraint program m, the algebraic quantum 
gravity program [T3j, the deparametrized models HMS1 in the canonical setting, the spinfoam 
program [2D] in the covariant framework and also some toy models |211124] in which an alternative 
quantization strategy of the Dirac algebra is applied. 

Concerning the technical difficulties, the Hamiltonian constraint is composed of two terms: 
the Euclidean part and the Lorentzian part. Both are non polynomial in the canonical variables, 
specially the second term that involves a double Poisson bracket of the Euclidean part with the 
volume and has a complicated form in terms of Ashtekar variables. A clever way to tame the non 
polynomial character of the constraints is using “Thiemann’s trick”, i.e. replacing the classical 
non polynomial functions by Poisson brackets of polynomial functions with the volume and of 
the Euclidean part with the volume. Once promoted to an operator the resulting expression 
comprises several commutators containing the volume operator |25fl27j . While this procedure 
helps to bypass the non polynomial character of the constraint, the resulting operator however 
is not self-adjoint and the explicit calculation of the Hamiltonian action is impossible because 
the volume operator present in the final expression has no explicit spectral decomposition. The 
partially formal result is already an extremely involved expression [[251 [29] . 

In this work we present another proposal for quantizing the Hamiltonian constraints. The 
first change is already in the classical formula for the scalar constraint. It is the sum of terms 
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proportional to the Euclidean scalar constraint and, respectively, the Ricci scalar of the three 
metric tensor [30]. Our aim is to implement the dynamics in the quantum model of gravity 
coupled to a free scalar field m- The construction is conceptually based on the recently in¬ 
troduced “intermediate” Hilbert space 77 v tx [3T] that is preserved by the obtained Hamiltonian 
operator, raising hope for a well defined evolution operator with satisfactory properties, e.g. self 
adjointness. 

The developed regularization is based on a concrete implementation of a proposal first ap¬ 
peared in [32] concerning the Euclidean constraint, and the use of the curvature operator intro¬ 
duced in [ 33] to deal with the Lorentzian part. The paper is organized as follows. In section 
[II] we review the classical model of gravity minimally coupled to a scalar field; in section III 


we review the loop quantum gravity construction, present the regularization of the Hamiltonian 
and discuss the quantum operator and its properties; then we close in section IV with some 
conclusions and outlooks to further developments of this program. 


II. CLASSICAL THEORY 

Considering gravity minimally coupled to a scalar field in the standard ADM formalism 3~il . 
the theory is set as a constrained system for the standard canonical variables q a b(x) and 4>(x), 
respectively the metric and the scalar field on a 3 d manifold E with conjugate momenta p ab (x ) 
and 7 t(x). The analysis shows that the vector constraints C a (x) and the scalar constraints C(x) 
in this model are expressed in terms of the vacuum gravity constraints, Cf r (x) and C gr (x), and 
the scalar field variables as follows 

C a {x ) = C| r (x) + 7T (x)(j) t a(x), (II.1) 

C(x) = C gT (x) + ]- 71 -M= + \(f b {x)^, a {x)^ !b [x)yjq{x) + V((/))>/q(x), (II.2) 

2 v q(x) 2 

where q is the determinant of the metric q a b- 

With the Ashtekar-Barbero variables (A l a ,Ef) (i = 1,2,3) used in LQG, 

{tf a {x),E]{y)} =SKf3G5 b a 8)5{x,y) (II.3) 

=0 = {E^x),E b {y)} (II.4) 

where G is Newton constant and (3 is the Immirzi parameter, additional constraints - the Gauss 
constraints generating Yang-Mills gauge transformations - are induced: 

G\x) =d a E“ + e ij k AiE a k . (H.5) 

The field A l a is identified with an su( 2 )- valued differential 1 -form 

A = A i a r i ®dx a (II. 6 ) 

while the field Ef with an su(2)* vector density 

E = (II.7) 

1 dx a 

where t\, T 2 ,73 Gsu( 2 ) is a basis of su( 2 ) such that 


= S t 


ij- 


— 2 Tr nTj 


(II.8) 
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A solution by points in the phase space (A z a , Ef, cj>, 7r) must satisfy all the constraints: 

G^x) = 0 Ca(x) = 0 C{x) = 0. (II.9) 

In terms of the Ashtekar-Barbero variables, the gravitational part of the scalar constraint reads 


c(x) = -w^g{ — 71 


(e ijk E<?{x)E b (x)F k ab {x) 


det E[x)\ 


+ (1 + /3 2 )y/| det E{x) \ R{x 


(II. 10) 


where R is the Ricci scalar of the metric tensor q ab on S related to the Ashtekar frame variable 
by 


q ab = 


EK 

| det E\ 


( 11 . 11 ) 


The first term of C gr is usually related to the Euclidean scalar constraint 

e ijk E?(x)E b (x)F* b (x) 


C 


Eucl 


1 


167 tG/3 2 det E(x)\ 


( 11 . 12 ) 


To construct a quantum theory mainly two strategies can be adopted. The first is to promote 
the whole set of constraints to operators defined in an appropriate Hilbert space and look for 
the states annihilated by the constraints operators to build a Physical Hilbert space. The 
second, that we consider in this work, is to deparametrize the theory classically then quantize. 


The deparametrization procedure starts with assuming that the constraints (II.9) are satisfied, 
hence we can solve the vector constraints for the gradient of the scalar field, 


Cf 


7r 


(11.13) 


and then use this condition in (II.2) to solve it for 7r: 


7T 2 = 


Vq (- (<? gr + vwm ± \[ (C^ + vwm 2 - . (11.14) 


In case of vanishing potential 


v(<f>) = o, 


(11.15) 


which is our assumption in the rest of this article, equation (11.14) represents the deparametriza¬ 


tion of the system with respect to the scalar field, which can be seen as an emergent time. Note 
that in this case, on the constraint surface, it is necessary to have 


C sr (x) < 0. 


(11.16) 


The sign ambiguity in (11.14) amounts to treating different regions of the phase space, namely 
for + and — respectively 


7T 2 > / < q ab (x)c/> t a(x)(/> }b (x)q(x). (H.17) 

We choose the phase space region corresponding to + and >. It contains spacially homo¬ 
geneous spacetimes useful in cosmology. Then, the scalar constraints can be rewritten in an 
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equivalent form as 


C'(x) = 7r(x) 4= h(x), (11.18) 

h := y -y/qOr + Jqsj (C^) 2 - q ab C^c[ r . (11.19) 

We will also restrict ourselves to the case of 

7r(x) > 0, (11.20) 

although technically there is no problem in admitting both signs in the quantum theory. 

The constraints C' commute strongly, 

{C'{x),C\y)} = 0, (11.21) 

implying |16| 

{h(x),h(y)} = 0. (11.22) 

In this case a Dirac observable G on the phase space would satisfy 

{G, G\x )} = {G, C a (x )} = {G, C\x)} = 0. (11.23) 

The vanishing of the first and second Poisson brackets induce gauge invariance and spatial 
diffeomorphism invariance respectively. The vanishing of the third Poisson bracket is equivalent 
to writing 


dG 
d <f> 


{^,7t(.t)} = {G,h(x)h 


(11.24) 


III. QUANTUM THEORY 
A. The general structure 

The quantization of gravity coupled to a massless scalar field was performed in mm- While 
the derivation was partially formal - the existence of the operators Uf 1 is assumed at some stage 
- the result is expressed in a derivable way by elements of the framework of LQG: 

• The physical Hilbert space 7i is the space of the quantum states of the matter free gravity 
which satisfy the quantum vector constraint and the quantum Gauss constraint. 

• The dynamics is defined by a Schrodinger like equation 

5* - -s** (nu > 

where t is a parameter of the transformations 

ip !->• ip + t . 


• The quantum Hamiltonian 


d?x\J—2yJ q(x)C sr (x ) 


H 


(III.2) 
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is a quantum operator corresponding to the classical observable 


H 


d 3 x\J—2y/q(x)C&(x). 


(III.3) 


This operator could be defined by using already known operators \Jq{x) and C& T (x), as outlined 
. However, the observable sjqC gr written in terms of the Ashtekar-Barbero variables reads 


m 


— y / |det£ , (x’)|C gr (x) = 


1 


16 irf3 2 G 




k E“{x)E b Ax)F* h {x) + (1 + /3 2 )\detE(x) \ R(x). (III.4) 


The de nominator y/\ det E(x)\ present in (II. 10) disappears in (III.4). Moreover, the formula 
(III.85) below for ^/|det-E(x)| R{x) expressed in terms of the quantizable observables (holonomies 
and fluxes) also contains the same denominator, which again disappears after using the formula 


(III.4). That coincidence of reductions motivates us to quantize the expression (III.4) for h(x 


directly. 


B. Kinematical Hilbert space 


The kinematical quantum states in LQG are cylindrical functions of the variable A, i.e., they 
depend on A only through finitely many parallel transports 


h e [A] 



(III.5) 


where e ranges over finite curves - we will also refer to them as edges — in S. That is a kinematical 
quantum state T has the form 


*[A] = 'iP(h ei [A},...,h en \A\) (III.6) 

with a function ip : 811(2)" — > C. Here, for every edge we choose an orientation to define the 
parallel transport h ei [A]. 

To calculate the scalar product between two cylindrical functions and \E f/ we find in E an 
embedded graph 7 " = {e'{. ..., e"„}, such that both functions can be written as 1 

T[A] = ip(h e AA],.. .,h e " 

1 n 1 

T'[A] = /(/t e "[A],..., V 

1 r 

The scalar product is 

(*!*') = ..-.*■)• (ni.8) 

We denote the space of all the cylindrical functions defined as above with a graph 7 by Cyl 7 
and, respectively, the space of all cylindrical functions by Cyl. The kinematical Hilbert space 
%kin is the completion 


„M). 


(IH-7) 


^kin = Cyl 


with respect to the Hilbert norm defined by (III. 8 ). 


(III.9) 


1 The existence of •y" is ensured by assuming the analyticity of E and of the edges of the graphs 0] . 
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Every cylindrical function / is also a quantum operator 

(f(A)*')[A} = f[A]*'[A]. (III. 10) 

A typical example is 

f(A) = D^\{h p (A)) (III.11) 

defined by a path p in E, a half-integer j = 0, 1, |,..., the corresponding representation 

D {j) : SU(2) -> Unitary(^ (i) ) (III.12) 

and some orthonormal basis vi, ...,V2j+i £ T~L^\ 

D (j)a b (g) := (v a \D^(g)v b ) nU) . (III.13) 

Note that a connection operator “A” itself is not defined. 


An operator J x r e i£, which is naturally defined in this framework, is assigned to a triple 
(x, £, [e]), where x £ E, £ £ su(2) and [e] is a maximal family of curves beginning at x such 
that each two curves overlap on a connected initial segment containing x. To define the action 
of J x [e\i on a function T £ Cyl, we represent this function on a graph such that ej £ [e]. The 
action is 

J x [e]^ = ^{h ei e^,h e2 ,...,h eri ). (HI-14) 

0,6 e=0 

For ( = T*, it is convenient to introduce a simpler notation 


J 


x,e,i 


:= J 


x\e\Ti■ 


(III.15) 


The field E?(x) is naturally quantized as 

= jimiEtix)} 


8tt/31p 5 

* SA i( x ) 


T[A], 


(III.16) 


Given an edge e : [foAi] ~► E, and a function / £ C(SU(2)), the variation is given by the 
following formula 


5 

5A i a {x) 


f(h e (A)) 



dte a (t)5(x,e(t)) {h e , tl ,t(A)T l h e ,t,to(A))^ 


d 

9 9b 


f (9^\g=h e (A) ’ 


(III.17) 


where by /i e ,t 1 ,t(A) (respectively, /i e ,t,t 0 (A)) we mean the parallel transport with respect to A 
along e from the point e(t) to e(t\) (e(fo) to e(t)), and by the partial derivatives with respect 
to group elements we mean 


f(ge^) ■= (gOi-Am- (hi.18) 

=o 9 9 b 

Smearing along 2-surfaces leads to well defined operators in T-L kin . Let S C E be an oriented, 
2 -dinrensional surface, and 


d_ 

de 


£ : S —> su(2) 


(III.19) 
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be a smearing function. The flux corresponding to E is 

Pst ■= [ Ux b Adx c e abc C(x)E?(x). (III.20) 

J s 

The quantum flux is a well defined operator 

Pst = 47 TGj2e(x)J2^s(e)J x ,e,i (III.21) 

xeS e 

where e runs through the classes of curves beginning at x. and 

«s(e) = -1,0,1, (III.22) 

depending on whether e goes down, along, or, respectively, up the surface S. A generalized 
function £ may also involve parallel transports depending on A. A typical example is 

£(x) = Ad(h p(x) (A))(, C € su(2) (III.23) 

where 

x i-)- p(x) (III.24) 

assigns to each point x a path p(x), h p ^ x )(A) is the parallel transport, and Ad is the adjoint 
action of SU(2) in the Lie algebra su(2) 

Ad( 5 )C = gtg- 1 . (III.25) 

In conclusion, the operators compatible with the LQG structure of "Hkin are (functions of 
the) parallel transports and fluxes. 

The quantum Gauss constraint operator reads 

Gi(x) = ^2 3( x , y) ^2 (III.26) 

y&Tj e at y 


Solutions in Cyl to the Gauss constraint 

Gi(x)tf = 0 (III.27) 


are functions such that 

/(A) = /{g^Ag + g^dg), for every j€ C 1 (E,SU(2)). (III.28) 

We denote their algebra, subalgebra of Cyl by Cyl G , and the corresponding subspace of "Hkin 
by 77^ n . A dense subspace of 77{(( n is spanned by the spin network functions. A spin network 
function is defined by a graph 7 = (ei,..., e n ), half integers (non zero) (j \, ...,j n ) assigned to the 
edges and intertwiners (ti,..., L rn ) assigned to the vertices (ui, 

T(A) = D^XiheM)) - D^X{h en {A)){A® ... ®Lmt\± n - ( IIL29 ) 

Each i a is an invariant of the tensor product of the representations assigned to the edges ej 
whose source is v Q and the representations dual to those assigned to the edges whose target is v a . 

Given a graph 7, we denote by Cyl^ the space spanned by all the spin network functions 



defined on this graph, and 


:= CylG (III.30) 

To define the orthogonal decomposition of the space of the Gauss constraint solutions we need 
to admit closed edges, that is edges for which the end point equals the beginning point, and closed 
edges without vertices (embeddings of a circle in £). In the case of an edge without vertices, 
we choose a beginning-end point arbitrarily in the definition of the spin network function. On 
the other hand we do not count those graphs that can be obtained from another graph by the 
splitting of an edge. Then the space of all the solutions to the Gauss constraint can be written 
as the orthogonal sum 


^kin = ©«? (III.31) 

7 

where 7 ranges over all the un-oriented graphs defined in this paragraph. 


C. The vertex Hilbert space 


Every analytic diffeomorphism / E DifP(E) defines a unitary operator Uf : %kin —> V- v in , 

U f V[A] = ^[f*A\. (III.32) 

Given a graph 7 consisting of edges and vertices 


Edge ( 7 ) := {ei,..., e n }, Vert ( 7 ) = {vi,...,v m }, 


the action of Uf on a cylindrical function (III. 6 ) reads 


Uf*[A] = iP(h f{ei) [A},...,h f{en) [A\), 


(III.33) 


where for the parallel transport along each edge /(ej) we choose the orientation induced by the 
map / and the orientation of e/ chosen in (III. 6 ). Smooth diffeomorphisms map analytic graphs 
into smooth graphs, therefore their action is not defined in our Hilbert space H \d„ ■ Suppose, 
however, that given a graph 7 , a smooth diffeomorphism / E Diff°°(E) maps 7 into an analytic 
graph. Then (III.32) and (III.33) define a unitary map 


U f 


n 


/( 7 )' 


(III.34) 


The idea of the vertex Hilbert space of m is to construct from elements of the Hilbert space 
n partial solutions to the vector constraints, by averaging the elements of each of the sub¬ 
spaces T~L 7 with respect to all the smooth diffeomorphisms Diff°°(E)y ert ( 7 ) which act trivially in 
the set of the vertices Vert( 7 ). Denote by TDiff°°(E ) 7 the subset of Diff°°(E) which consists of 
all the diffeomorphisms / such that /( 7 ) = 7 and Uf acts trivially in 'H~j . 

Denote by Diff“(E)v ert ( 7 ) the set of those elements of Diff°° (£)y ert ( 7 ) which preserve the 
analyticity of 7 . The maps % 7 —> obtained by the diffeomorphisms in Diff“(E)y ert ( 7 ) are 

in one to one correspondence with the elements of the quotient 


D 7 := Diff“(E) Vert(7) /TDiff(E) 7 . 


(III.35) 


Since D 7 is a non-compact set and we do not know any probability measure on it, we define the 
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averaging in Cyl*, the algebraic dual to Cyl. Given T G % G , we turn it into (*I'| G Cyl*, 

<*| : VGCyl G :=®Cyl? ^ (tf|tf'), 

7' 


and average in Cyl*, 

^ E ( IIL36 ) 

7 [/]er> 7 

where 


iV 7 = |Sym 7 | 

Sym 7 := {/G Diff”(S) {ll) ,.., m} : /( 7 ) = 7 } / TDiff(E) 7 . (III.37) 

The resulting 77 ('L) is a well defined linear functional 

7 ?(T) : Cyl G -> C 

because given \E r/ G Cyl G , only a finite set of terms in the sum contribute to the number r/(\I')(T / ). 
Hence we have defined a map 


% G 9 T i-7 r/O) G Cyl* 


for every embedded graph 7 . We extend it by linearity to the algebraic orthogonal sum 

V ■ ^ki„ —> Cyl*. 

The vertex Hilbert space 'H G X is defined as the completion 


ng x := 7 (Cyin^ G n ) 

under the norm induced by the natural scalar product 


pIPI) 

(III.38) 


(III.39) 


(r / (^)|7 7 (V)) := 77(®)(®0- 


(III.40) 


It has an orthogonal decomposition that is reminiscent of (III.31): Let FS(£) be the set of finite 
subsets of £. Then 


n_/G 

^vtx 

0 ^ 

£ 

© 

II 

(III.41) 


vsfs(s) 


njG 

tty 

0 

11 

(III.42) 


[7]e[7(V)I 


cG 

:= 

(III.43) 


where 7 (V) is the set of graphs 7 with vertex set V = Vert( 7 ), [ 7 (H)] is the set of the Diff(E)y- 
equivalence classes [ 7 ] of the graphs 7 G 7 (H) and 5 G is the subspace 5 G C Lf G of the elements 
invariant with respect to the symmetry group Sym 7 . Importantly, 

V = S G ^ S G (III.44) 

is an isometry. The orthogonal complement of 5 G in % G , on the other hand, is annihilated by 
V- 
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The Hilbert space 
It is defined by 


carries a natural action of Diff UJ (S), which we will also denote by U. 


U f rj(V) :=ri(U f V), f G DifF(£) (III.45) 

A short calculation shows that Uf is unitary and maps 'Ey to 'H ( ^ V] in the decomposition 

( jlll.41| ). 

Each subspace <S£ consists of DifP(£)y ert ( 7 ) invariant elements. In this sense, they are partial 
solutions to the quantum vector constraint. They can be turned into full solutions of the quantum 
vector constraint by a similar averaging with respect to the remaining Diff(£)/Diff(£)y ert ( 7 ) [31]. 
We denote the space of those solutions 

D. The Hamiltonian operator 

In the Hilbert space 'H^. x we will introduce (derive) an operator 

H = j d 3 x^J-2^(x)C&{x), (III.46) 

where we have 

- 2 \J | detE(x) |C gr (x) = (e ijk E“(x)E b j (x)F^ b (x) + (1 + f3 2 )\detE(x)\ ®R(x)?\ (III.47) 

In order to dehne the corresponding operator, we need to consider how to regularize and quantize 
an expression of the form 


J d 3 x f (x) y/a 2 (x) + b 2 (x) (III.48) 

where / is a smearing function defined on £ while a(x) and b(x) are functionals of the fields A\ 
and Ef. 

Introducing a decomposition of the manifold £ into cells A, the integral can be approximated 
as 


d 3 x /(x)y / Q 2 (x) + b 2 (x) = Jd 3 xf(x)a(x) S j + (^j d 3 xf(x)b(x) S j + 0{e a), 


(III.49) 

where for every cell A, denotes the coordinate volume of A. If the integrals J A d 3 x f(x)a(x) 


and f A d 3 xf(x)b(x) can be quantized as well-defined operators, equation (III.49) then shows 
how to dehne the operator corresponding to J s d 3 x f(x)^a 2 (x) + b 2 (x). Equation (III.49) is 
the basis of our construction of the operator ( III.46| ). 

In our case, the operators corresponding to a(x) and b(x) themselves will be available, and 
will have the general form 


a[x = 


£*( 


x , v )a v 


v&Yj 


b( x ) = Tj(x,v)by, 


(III.50) 


where the operators a v and b v , when applied to a spin network state dehned on a graph, have 
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a non-zero action only if v is one of the vertices of the graph. In this case, the operator 


d 3 x /(x) i/a 2 (x) + b 2 (x) 


(III.51) 


can be defined simply by inserting a(x) and b(x) into the right-hand side of equation (III.49). In 
this way one obtains an operator, whose restriction to the space of spin network states defined 
on a given graph 7 takes the form 


d 3 x f{x) \Ja 2 (x) + b 2 (x) 


= ^f{v)\Jd 2 + b 2 . 


(III.52) 


In other words, our regularization gives 

\]a 2 (x) + b 2 (x) = E 5{x,v)\ja 2 + b 2 . 


(III.53) 


v&T, 


1. Euclidean part 


We start with the quantization of the Euclidean part of our Hamiltonian (see (11.12)). In 

of a 

H Eac \x ) : = 


equation (III.49), the role of a(x) is now played by the function 

1 


,^/e ijk E9(x)E^x)F^(x). 


^8irGp 2 

Consequently, we consider the quantization of the integral 

[ d 3 xf(x)H Eucl (x ) 


(III.54) 


(III.55) 


(where an arbitrary smearing function / has been introduced). 

According to the general framework of LQG, we need to express the integral in terms of 
parallel transports h e and fluxes Ps ,i- The easiest example is to consider the Riemann sum for 
this integral obtained by considering a cubic partition V e of E into cells □ of coordinate volume 
e 3 

e 3 E ^e ijk EUxn)E!>(x n )FZ b (x n ), (III.56) 

□ 

and to distribute the e suitably 

E -^ Xn ) \J e b' fc ( e2E i (*□))( e2 Ej (x n )) {e 2 F^ b (x n )). (III.57) 

□ 

For each cube □ denote by xu the center, by a = 1, 2, 3, three sides x a =const (for each a 
there are two, choose any one and orient such that the following is true). Moreover, for every 
x £ □, denote by pn(x) the line from xn to x G □. Then, we have 

e 2 £“ = Ps h ,i + o{e 2 ) 
e2 Eab = tabcihgg^ )^ + o(e 2 ), 


(III.58) 
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where by Ps a ,i we mean Pg^ of (III.20) with 


S = Sg, and £(x) := h pn(x) Tih p ^ x) , (HI-59) 

h pn ^ standing for the parallel transport (with respect to a given field A) along pu, and for a 
SU (2) element h we define 2 


(/j fc )(0 = Tr 

w? 



(III.60) 


with W] ; = iy/l(l + 1)(2 1 + 1). In this way we write the original expression in terms of fluxes 
and parallel transports (as a limit), 


J d 3 xfH EucYPe 


1 

tG/3 2 


y! f ( x °) \J e ijkPs^iPs^ J e “fcc(^95g)^\ 
□ 


(III.61) 


in the sense that 

J d 3 xfH Enclv * -> J d 3 xfH Euc \ (III.62) 

in the limit e —> 0 when we refine the partition (□ —> •). 


More generally, we regularize the integral by using a partition V e which consists of: 

• an e-dependent cellular decomposition C e of X; 

• assigned to each cell A € C £ : 

— a point xa inside A; 

— a family of 2-surfaces C dA, 1 = 1, ...,71a; 

— a family of paths pa{x) labeled by points x G dA, each going from xa to x; 

— a family of loops a A , K = 1,..., 777a; 

— a family of coefficients ka ijk\ 


such that the following functional 


[d\fH E ^-(A,E) := ‘ X fM.Ak £ ^AlcPy-psZMAfi (III.63) 

J V 87r(j P A V IJK 

approaches the Euclidean Hamiltonian, 


d 3 xfH EucYPt {A, E) 


e —>0 


d 3 xfH 


Eucl 


As in the cubic example, by P s i^ i we mean Pg^ of (III.20) with 


S = S A , and f(x) := h p ^ x) Tih p ^ {x) 


(III.64) 


(III.65) 


Equation (|III.60[) is obtained using the relations 


D m {h) = h°l m +h k (r k ) m , Tr(ri) (0 =0, Tr {{ri) (l) (r k ) (l) ) = ^-6*. 
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Each term £ijk K ATJKPs' A iPs£j ^ l a K )^ gi yes Hse to a well defined operator in 77km 

tijkKAUK : P s i ,iP S i,j () W: (HI.66) 


which depends on the ordering of the operators, symbolized by : 

In this way we obtain an operator 

J d 3 x/(x)77 Eucl ^(z) := ^/ 8 ^ G/32 J2 f(. x A)\Jei jk HAiJK-P s i A ,iP S i j(h*jr) W = (III.67) 

which depends on the partition V e and is well defined in 77km- However, as we refine the partition 
V e . the operator family does not converge to any operator in 77km- This is a well known problem 
in LQG and it does not have a solution in the kinematical Hilbert space 77 km- 

A way out is to consider the dual action of the regulated operators / d 3 xf(x)H Encl ^ e (x) in 
the Hilbert space 77(^ x . That was done for the (formally regularized) operator C gr in [31]. As 
it is explained therein, and those arguments apply also in the case at hand, a limit as e —)• 0 
exists upon several conditions about the partitions V e . To begin with, we adjust the partitions 
individually to each subspace 77 7 in the decomposition (III.31). Secondly, a successful partition 


has to have a suitable diffeomorphism covariance in the dependence of the partitions on 7 and 
on e. 

The outstanding problem though, is the dependence of the result on choices made. There are 
many partitions which satisfy the conditions. The resulting operator carries a memory of the 
choice of V e , for example on the adjustment of the fluxes to graphs. To restrict that ambiguity, 
we study first the straightforward quantization of E b . 


Let T G Cyl be as in (III. 6 ). Assuming E. 


x = T 


i 5A l a (x ) ’ 


we obtain 


e ijkF a b( x )\ ■ 


h 


: ^(h ei (A),...,h en (A)) = 


ij 5A i a {x) 5A{(x) 

Y f [ dtdt'S(x,e I (t))6(x,e r (t l ))F^ b (x)ei(t)e b I ,(t ^ 


e/,e 7 / 


/to Jt' 


A f \ ^ d d 

{he I ,t 1 ,t( A ) T ^e 7 ,t,t 0 (^)) B [heph ep>t ',t' 0 {A)) B , '" ,gn ^\9i=h ei (A), gi ,=h ei ,(A) 

(III. 68 ) 


Certainly the product of the two Dirac delta distributions is ill defined at some points x, ei{t) and 
eji (t'). However, we can precisely indicate those points at which the expression is identically zero. 
To begin with, the product S(x,ei(t))5(x,ep(t')) vanishes except for the triples (x, e/(7), ep(t’)) 
such that 


x = ej(t) = ep{t'). (III.69) 

Secondly, the factor F^ h (x)e 0 f(t)e b T ,(t r ) is not zero only if 

e/C t) Merit'). (III.70) 

Those two conditions are satisfied simultaneously only if x coincides with one of the vertices 
v of 7 and the edges e/ and ej> intersect transversally at v. Suppose that x = v £ Verify), the 
edges e/ : [7 q,7i] S and ep : \t' Q . t\\ -A S intersect transversally at v, both oriented to be 
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outgoing. Then the corresponding contribution comes only from 

x = v, tj = to, tp = t' 0 , 


(III.71) 


and it is 


(KM)t‘) a b {he,, (Ay) 


A' d d 


B' 


A' 


dgiB dgp B 


-,9n) 


gi=h ei (A),gp =h ep {A) 

(III.72) 

modulo the ill defined factor (5(v,v)) 2 which has to be regularized. Our regularization is also 
expected to replace F^eiey by a parallel transport h en , along a loop ew assigned to the two 
(segments of) edges. Finally, diffeomorphism invariance implies that each vertex v and a pair 
of transversally intersecting edges ej and ey at v contribute the same operator as any other 
diffeomorphism equivalent triple v ', e) and e),. 

We are now in a position to formulate assumptions about the construction of the partitions 
V e adapted to a graph 7 , as shown in ED, and the assumptions about the assignment of the loop 


a A used to regularize the connexion curvature Ff ,, in order to guarantee the diffeomorphism 
covariance of the final operator. 

Given a graph 7 = (ei, ...,e n ) of Vert( 7 ) = (tq, ... ,v m ), in order to spell out the conditions 
it is convenient to split each edge into two segments and orient the new edges to be outgoing 
from the vertices of the original graph 7 . Denote the resulting graph by 7 ' = (e^, e' 2n ) and 
its vertex set Verify') = (iq,..., v m , v ' m+l ,..., v' m+n ). The assumptions are as follows: 

Requirement 1. 

• each cell A contains at most one vertex of the graph y'; 

• if v £Vert( 7 ) and v £ A , then 


— xa = v; 

— to each edge e\ there is assigned a surface S A C dA intersecting the edge transver¬ 
sally (there may he surfaces in dA not intersecting any edge); 

— to each ordered pair of edges e 1 and ef meeting transversally at v there is assigned a 
loop a A oriented according to the order of the pair (e), e'j ). Hence we denote kaijk 
as kaijij; 

— for edges e) and e'j of 7' meeting transversally at v k-aijij is not zero; 

— for edges e) and e'j of 7' meeting tangentially at v At a ijij = 0; 

— for edges e) and e'j of 7' meeting transversally at v the corresponding loop ajj, in 
the limit e —> 0, is shrank to v in a diffeomorphism invariant way; 

— if A does not contain an edge of 7 but it contains a segment of an edge then, by 
splitting the edge and reorienting its segments suitably, we turn that case into the 
case of A containing a 2 -valent vertex; 

• the value of non vanishing At a ijij is an overall constant Ati(u) depending on the valence 
of the vertex but independent of A, I, J. 

Concerning the prescription for the assignment of the loops ajj - we call them special loops - 
which are created by the Euclidean part of our Hamiltonian operator, we wish the construction 
to satisfy the following requirements: 

- The loop added by the Hamiltonian should be attached to the graph according to a 
diffeomorphism invariant prescription El ED- This property allows the operator to be 
well defined on the space H^ tx . 
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- It should be possible to distinguish between loops attached to the same vertex but asso¬ 
ciated to different pair of edges, and between loops attached to the same pair of edges 
by successive actions of the Hamiltonian. This property makes it possible to define the 
adjoint operator on a dense domain in TL^ tx , and consequently to construct a symmetric 
Hamiltonian operator. 

Consider a vertex v of the graph 7 defined above and a set of links {e/} incident at v. In 
order to satisfy the first requirement, we use a construction that was introduced in [ 36] and was 
presented in a work of T. Thiemann [9]. The construction consists of two parts. Firstly, to each 
pair of links ej and e j incident at v, we define an adapted frame in a small enough neighborhood 
of v. Then we require that the loop aqj, associated to the pair (e/,ej), lies in the coordinate 
plane spanned by the edges ej and ej. The choice of the adapted frame is based on the following 
lemma: 

Let e and e' be two distinct analytic curves intersecting only at their starting point 
v. Then there exist parameterizations of these curves, a number 5 > 0, and an 
analytic diffeomorphism such that, in the corresponding frame, the curves are given 

by ' 

(a) eft) = (i, 0, 0), e'(t) = (0, t, 0), t G [0,5] if their tangents are linearly indepen¬ 
dent at v, 

(b) eft) = (t,0,0), e'(t) = (t,t n , 0), t G [0,5] for some n > 2 if their tangents are 
co-linear at v. 

We will call the associated frame a frame adapted to e, e'. 

To carry out the second part of the construction, we need a diffeomorphism invariant pre¬ 
scription of the topology of the routing of the loop ajj. In other words, the plane in which the 
loop lies should be chosen in a way which is diffeomorphism invariant, and which does not cause 
the loop to intersect the graph 7 at any point different from the vertex v. The choice that ajj 
lies in a small enough neighborhood of v guarantees that the loop cannot intersect any edge of 
7 except the edges incident at the vertex v. Then the routing of the loop in that neighborhood 
is achieved through the prescription given in [9] (and which we do not repeat here). 

Now let us turn to the second requirement, which is crucial in order to have the possibility of 
defining a dense adjoint operator that allows one to construct symmetric Hamiltonian operators, 
and eventually to provide self-adjoint extensions. To state the prescription that satisfies the 
second requirement, we need to define the order of tangentiality of an edge at the node. This is 
defined as follows. Considering the vertex v and the edge ej. we denote by kjj > 0 the order of 
tangentiality of e/ with another edge ej incident at v. If the edges ej and ej are not tangent at 
v, we understand that k jj = 0. The order of tangentiality kj of the edge ej at the vertex v 

kj = max kjj (III.73) 

ej at v 

i.e. as the highest order of tangentiality of the edge ej with the remaining edges incident at v. 

The element which completes the prescription of the special loop according to the two re¬ 
quirements is now stated as follows: 

Requirement 2. 

The special loop au is tangent to the two edges ej and ej at the vertex v up to 
orders kj + 1 and kj + 1 respectively, where kj(> 0) and kjf> 0) are respectively 
the orders of tangentiality of ej and ej at the node. 
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This property indeed makes a loop attached by the Hamiltonian to a given pair of edges 
perfectly distinguishable from any other loop at the same node. 

To summarize, the prescription for assigning a special loop to a pair of links incident at a 
vertex is to choose the loop to lie in the coordinate plane defined by the frame adapted to the 
pair of edges, then to follow a specific and well defined routing of the loop described in [9J, and 
finally to impose the tangentiality conditions introduced above. With this prescription, the 
loop assigned to a pair of edges is unique up to diffeomorphisms. 


In consequence, given a graph 7 and the auxiliary graph j 1 obtained by the splitting, the 
contribution from a cell A containing a vertex v reads 


K i( v )^2 € ijk e(e'j,ej) :J, 


J ?! 




i(h k u)) 0) ' : =■ 


(e) 


(III.74) 


i,j 


where e (e'j, e'j) is 0 if e) and e'j are linearly dependent or 1 otherwise. This operator maps 


--(e) 

m ■■ Cyi 7 


Cyly// 

l" = 7 


(III.75) 


IJ 


Considering all the graphs we combine the operators into a single e-dependent operator 

-M 


H* : Cyl ->■ Cyl . 


(III.76) 


In order for the operator (III.76) to be cylindrically consistent, we should have k\(v) = k,\, 
an overall constant independent of the valence of the vertex v. 3 However, since our goal at the 
end is to implement this operator in the gauge invariant Hilbert space, we can equally well define 
the operator by proceeding with the regularization directly on the spaces Cyl G orthogonal to 
each other. In that case the question of cylindrical consistency does not arise, and we may allow 
the possibility that ki(v) depends on the valence of the vertex. 


In this way we have determined the action of an operator (III.55) up to a value of k±(v) 
(constant or not), assuming the conditions (111) and (12|) . This operator passes naturally to Cyl G 


H.} {€) : Cyl G Cyl 


G 


(III.77) 


As we refine the partition by e —> 0, the loops ajj are shrank to v. However, the e-dependent 
operator ( ] defined by the duality * in 77 G X (on a domain that includes r/(Cyl G )), 


7(4/) H7 ( ) 7/(tt) 


HZ (e) ) ? ? (T) : ^ 1 y ri(*)(HE*'). 


(III.78) 


is insensitive to the shrinking, as long as each loop ajj is shrank within the diffeomorphism 
class of 7 U ajj. Hence we drop the e label in the dual operator. It follows that the Euclidean 


On the other hand, every value of the constant Ki can be achieved by a suitable choice of the shape and size of 
the loops aij. 


3 
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part of the Hamiltonian is defined as 


f d^xf(x)H^\x) := — L=Y,f(. v )VHf ( IIL79 ) 

s V87TG/3 2 t&i 


In order to define the square root in this equation, one could choose a symmetric ordering of 
H E V . However, a symmetric ordering of the Euclidean term is not necessary for constructing the 
complete Hamiltonian, for which instead the square root of the sum of Euclidean and Lorentzian 
terms needs to be defined. 


2. Lorentzian part 


Following the strategy of quantization indicated by equation (III.49), we now introduce a 


second operator corresponding to the integral of the term \fqR again smeared with an arbitrary 
function / 


(i 3 xf{x)y/qR(x ) 


The construction of the operator is in two parts: first we write an approximate expression of the 
classical integral by implementing a cellular decomposition of the 3d manifold, characterized 
by a regulator e. Secondly, the regularized expression is promoted to an operator, which after 
taking the regulator limit, leads to a background independent operator acting in the Hilbert 
space of gauge invariant states. 

The aim is to construct an operator corresponding to the following function on the classical 
phase space: 


d 3 x /(x)-\/|det[E]| R = f d 3 x f{x)\J \/|det[E]| ^/jdetfEJIi?, (III.80) 


Consider a cellular decomposition c € e of the manifold E. The size of the cells is assumed to 
be controlled by the regulator e, in such a way that the coordinate size £a' of each cell A 7 E 


satisfies ca' < e. We can then write the integral (III.80) as a limit of a Riemannian sum over 
the cells A', 

j d 3 x f(x)\Jy/\det[E]\ y / |det[.E]|.R = lim ^ f{xA')J(^j d 3 x ^/\det[E}\^j (^j d 3 x ^\det[E]\R), 

(III.81) 

where on the right-hand side ia' denotes any point inside A 7 . 

Next we decompose each cell A 7 into ca' closed cells A, where a cell A has a boundary 


formed by a number ha of 2-surfaces (faces). In equation (III.81), we then approximate the 
integral of -y/|det[£ , ]| by a Riemannian sum over the cells A, and the integral of y^detDEJIi? by 
a regularized Regge action for an appropriate A-decomposition of A 7 , obtaining 


( d3 x f{x)yj \/|det[E]| y]det[E]|R = lini ^ /(®A')* ( Vqa(E) \ ( ^ Ra{E) J. 

: A'etf* \ \AcA' J Vac A' / 

(III.82) 
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The functionals q/\(E) [25] and Ra(E) are defined on the classical phase space as 4 

«o(A) 


Qa(E) := 


3! 


I^J^K 


e ijkP S i,i P Si,j P SK,h 


(III.84) 


Ra(E) := £ 

uCdA 


^ kP si >u J P si u ,k lw p si u ,j' p sl K ,k' 


( 2ir 

x-7r + arccos 

\ Oi u 


Vqa(e) 

5klP Si,„k P Si,„l 


(III.85) 


21 P, 


Ql 


IP, 


V 


where we use the following notation: 


- given A, the index 1 = 1, ...,tia labels the surfaces (faces) S A forming the boundary <9A 
of the cell A and u labels the hinges on that boundary (the 1-skeleton of the cell); 

- the symbols Si and Si stand for the two surfaces in dA that intersect at u; 


- the symbol P s i^ i represents the flux of the held Ef across S A , dehned in (III.20) with 

S = S ! a , and f(x) := h pA(x) Tih~^ x) 


and 


\] 5kk ' P sku* P s^ u ,k'i 


(III.86) 


(III.87) 


- k o(A) is a regularization constant depending on the shape of the cell A; 

- finally a u is a fixed integer parameter corresponding to the number of cells sharing the 
hinge u in the cellular decomposition c if e . 

Considering the coordinate size < e of the cell A, dehned such that the limit e —> 0 is 
equivalent to 6 a 0 , the functional q&{E) is such that v^qA^E) approximate the function 

|det[P]| at any point within the cell A, Va oc e A being the coordinate volume of A. Also, each 


term in the sum dehning Ra(E) (III.85), rescaled by L u oc ca that is the coordinate length of 
the edge u on the boundary of A, 


1 

Lu\ 


hijk p si j P Si k h'j'V P SL tU j' P SL u ,k' ( 2vr 

da' - . - - . -- 7 r T arccos 

V QA(E) v Qa( p ) V a u 


dklP Sj !U ,k P Sj u ,l 
21Pc,/ \\P qJ | 

1 °A.u" °A,u' J 


(III.88) 


approximate the function L U (E)Q U (E) in the limit €a —>• 0, where L U (E) and 0 n (P) are respec¬ 
tively the length of the hinge u and the dihedral angle at u in A expressed in terms of densitized 
triads. 


4 The functional qa(P) can be defined in a different way: 

«o(A) 


9a(E) := ]T 


IJK 


3! 


e IJK£ijkPsii P Si,j P S%,k 


(III.83) 


This definition would lead to a volume operator that is sensitive to the differential structure at the nodes, see 
[261 [27l|37]. 
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The sum over the cells A of the functional Ra(E) corresponds to the regularized Regge 
action [38] in 3d on A 7 , which is by itself an approximation of the function f A , d 3 x V\det[E}\R. 
We direct the reader to [33] for more details about the concepts of this construction. 


To continue the calculation from equation ( III.82| ), we assume that the cells A are chosen 
such that we obtain the same contributions q/\(E) and Ra(E) from each cell A, up to higher 
order corrections in e^' (equivalently, up to higher order corrections in e). Hence each sum over 
the cells A becomes the number of cells ca' times the contribution of the cell A, chosen as the 
cell containing the point xa> at which the smearing function / is evaluated. In this way we 
obtain 

J d 3 x f(x)yj \/|det[E]| ^et[E}\R = lim ^ f(x AcA ') C A '\j\/ qa( E ) r a( E )- ( IIL89 ) 


Let us now introduce the approximation 

f( x AcA') = T“ J2 /( x a), 

p- A / 


(III.90) 


AcA' 


which is an averaging of the values of the function / inside the cell A 7 , and which can be seen 
as a better approximation of the value of the function / inside the cell A 7 , in the sense that we 
are probing the function / in several points inside the cell instead of one point xa>- Inserting 


equation (III.90) in equation (III.89), we come to the result we are looking for: 


d 3 x f(x)\J v^det^l ^\det[E]\R = lim ^ ^ f {x a) \[\/qA{E) Ra(E) 

A'etft AcA' 

= lim ]T f{xA)\/y/q^E) R A (E), (III.91) 


Ac*f e 


where last step is achieved by combining the two sums over A 7 and A. 


Notice that the expression of Ra{E) in (III.85) contains an overall factor of (-y/ qA(E )) 


v-l 


This leads to a crucial simplification in the expression of yj^ jaJE ) Ra(E), namely, the factors 
of a/ qA{E) are canceled: 


Vqa(E) Ra(E) = J2 \ J 6 ii' e ijk p si u ,j P Si u ,k ti'j'k'Psi u ,f P Si u ,k' 


uCdA 


2vr 


J A U A ,ii’ J u A,u'>J u A,u 

5klP S^k P SiJ 


(III.92) 


XI-7r + arccos 

\ &IL 


J A ,u ’ 


2\P, 


cl I I PCj I 
A,” ^ A,u J 


In the quantum theory, this simplification implies that the volume operator will be absent from 
the Lorentzian part, and consequently from the whole Hamiltonian operator. The absence 
of the volume is an important technical advantage in the calculation of the action of the 
Hamiltonian. 

Before promoting this expression to an operator, we study the term €ij k P s i jPgj k ap¬ 


pearing in equation (III.92). This term approximates the classical function 


1 


£ijk£abcEjE k il (s) — ^hm^ ^4 e ijkPs^ u ,j P S^ u ,ki 


(III.93) 
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where s is parameterizing the curve u. 


Considering T in Cyl with a graph 7 = (ei,...,e n ), the straightforward quantization of 
(-Hk^-abcE^Euifis) by replacing E “ with j—f —- induces the factor 

J J lx J % 

e abc e?(t)e 6 J (t') (III.94) 

in the formal action of the operator on 4'. Similarly to the case of the Euclidean part (see section 


III D 1), this factor vanishes unless 


e,(t)tfej(0- 


(III.95) 


which means that the edges e/(f) and ej{t') are different (/ 7 ^ J) and transversal at their 
intersection point. In order to pass this property to the quantum operator, we introduce the 
coefficient k' aij , defined in the following, in the expression of \Jq/\(E ) Ra(E) and we write 


Vqa{E) Ra(E) := k 'a ij 2\J 5ii ' eijkP Si,u’j P Si, u ,k €i 'f k ' P Sk u ,j' P Si 

6klP sku* P s£, u ,i 


,k' 


(III.96) 


uCdA 


2ir 

-7T + arccos 

O',/. 


21 P, 


'si \\Pn 


J A ,u 


A ,u 


In order to promote the expression in (III. 96 ) to a quantum operator, we first need to set 
some requirements on the decomposition so that we adapt it to the functions in Cyl: given a 
T in Cyl with a graph 7 = (e\, ..., e n ) of Vert(7) = (ui, the requirements are as follows: 

• each cell A contains at most one vertex of the graph 7; 

• each 2-cell (face) on the boundary of a cell A, containing a vertex of 7 , is punctured 
exactly by one edge of the graph 7. The intersection is transversal and belongs to the 
interior of the edge; 

• if v <EVert( 7 ) and v £ A , then 

- xa = v, 

- s' AIJ is not zero only for edges ej and ej of 7 meeting transversally at v; 

• if A does not contain an edge of 7 but it contains a segment of an edge then, by splitting the 
edge and reorienting its segments suitably, we turn that case into the case of A containing 
a 2-valent vertex, 


A result that follows from the derivation of the curvature operator in [33] is that the value 
of non vanishing k' A jj is an overall constant K2{v) depending on the valence of the vertex but 
independent of A, I, J. This property is obtained from an averaging procedure used in order to 
remove the dependence on < ^ ?e . 

Having the quantum operators corresponding to P s i^ j s’, we are now able to define the 
quantum operator corresponding to q\{E) Ra(E) as 


Qa(E) Ra(E) ^ k! aij ~ y^u'^kPsi u ,k ^A'k'Psf 

u<ZdA 




(III.97) 


2 vr 

-7r T arccos 


hiPf 


ql 1 

°A,u^ 


2 5 kk ' P Si u ,k P Sik' 5 ll' P SiA P Sil' 
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Considering a cylindrical function \k 7 in the Hilbert space Cyl 7 , thanks to regularization 
detailed above we have 


Hi * 7 = 


{ Hi if A contains a vertex v of 7 ; 

0 if A does not contain a vertex v of 7 . 


where 


Hi: = K 2 (u)^e(e/,ej) i^ e/y 
i,j 

= - 2 - 

7,J 


ej 


2 vr 


&iV i,^ijkJv,ej,jJv,ej,k)^i'j'k' ^v,ei,j' Jv,ej,k ') 

^klJv,ej,kJv,ej,l 


— 7 r + arccos 


a 


v,e I: ej 


(III.98) 


(III.99) 


$kk' Jv,ei,kJv,ei,k' V I ^IV Jv ,e j ,lJv ,e j J 


with a t , je/iej an integer parameter, k 2 (u) is the averaging coefficient that depends only on the 
valence of the vertex v, and e (ey, ej) is 0 if ej and ej are linearly dependent or 1 otherwise. 


Now we can define an operator acting in Cyl G 



(fix f(x)y/\det[E]\ 


R := 





(III.100) 


that corresponds to the Lorentzian part of the (smeared) Hamiltonian. Consequently we intro¬ 
duce the operator Hi defined by duality * on ^vtx 


(1 + /3 2 ) f 
87 t^g y s 


(fix /(.'c)v / |det[£’]| R := 



(III.101) 


Notice that there is no ordering ambiguity in the operator Hi and therefore no ordering ambi¬ 
guity in Hi . 


3. The Hamiltonian operator & its properties 


At this level two operators H f' and Hi have been implemented on H v t x . Now we introduce 
the operator Hi defined as 


iijHi +(1 + /3 2 )Hi 

= J2e(e,e')[(Hi e ^y + (1 +/ 3 2 ) 


(III.102) 


e,e' 


where e, e! run through the set of germs of bounded one-dimensional sub-manifolds of £ incident 
to a point x, and we define Hi e e , as 


f 


•— 77i (xjCijk {h a ,)^ Jx,e,iJx, 


e ,3 ! 


(III.103) 


by choosing the simplest ordering for the Euclidean part. 
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a. A symmetric Hamiltonian operator 


The Hamiltonian operator represents the quantization of the classical Hamiltonian of the de- 
parametrized theory of GR coupled to a free scalar field. This final operator is required to 
be self-adjoint on some non-trivial domain in order for it to generate unitary evolution of the 
quantum system and for its spectra to admit a physical interpretation. Therefore, a first step 
toward achieving self-adjointness 5 of the Hamiltonian is to construct a symmetric operator. 

A symmetric operator H could be introduced as a combination of the operator Hj and its 

adjoint operator Hj . The later is defined on a domain f H J ] D S G , such that for every 


two states |\1/) and |^ ,/ ) in the space S G we have 




<*'| HJ \*) = (*\HZ\V) 


(III.104) 


where the bar stands for complex conjugate. 

For instance we choose the following definition for H 


H : 


H : = 


m t ) n 0 ( H? ) D S 


?G 


, 1 Wffj+gf. 
v / l6^^ V 


"Rvtx 


(III.105) 


b. Gauge and diffeomorphism invariance 

The operators h af , are J x ,e,i are both gauge invariant. Hence the Hamiltonian operator H is 
gauge invariant. 

Considering the group of smooth diffeomorphisms, the operator H is also diffeomorphism 
invariant, thanks to the regularization adopted and the averaging procedures involved in defining 
the curvature operator [33]. As a consequence of its gauge and diffeomorphism invariance, H 
preserves the space H G tx . 


c. Action of the symmetric Hamiltonian operator 


Looking into the action of H on a spin network state |y, {_)}, {/-}), one would like to express the 
matrix elements of this operator in terms of the quantum numbers labeling the states, namely 
the spins j and the intertwiners l. However in order to obtain such an expression, it is necessary 
to diagonalize the operator Hj under the square root in the definition of H. So far this aim 
of diagonalizing this operator has not been realized. Nevertheless, what can be achieved is the 
explicit computation of the action of Hff on a spin network basis st ate 1 7 , {. 7 }, {<■} ). We n eed 
first to calculate the action of Hff + (1 + fi 2 )Hf then from (III.78), (III.101) and (III.104) we 


can deduce the action of the dual operators, hence the act ion of H f and its adjoint. 

Moreover, having the state I 7 , {)}, {t}), we know from (III.102) that the action of H reduces 
to a sum of contributions, each coming from a pair of edges incident to the same vertex v of the 

graph 7 ; we call such a pair a wedge. The operators Hff ee ,, each associated to such a wedge in 


The self-adjointness of the Hamiltonian operator presented in equation (III.1051 is still an open question. 
However, we are able to construct self-adjoint extensions using other symmetrizations. 
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7 , can be introduced as 


HT. 


| 7 ,{.;},{,}) :=e(e,e') + (1 + /? 2 ) (i^ e>e ,)*] |7 ,{j}>W>- 


(III.106) 


Consequently we just need to compute the action of H^ e e , and H^ e e , leading to the total 

action of H^ ee , on I 7 , {j}, {/,}). The result is presented in the following. The calculations were 
done using the so-called graphical calculus, a framework that is briefly outlined in Appendix [A] 
along with the chosen notations and conventions. 

Consider again the spin network state I 7 , {j},{/.}} containing the n-valent vertex v (n > 1), 
from which originate the two edges e and e!. The actions of the Euclidean and Lorentzian parts 


of HT on the state 


v j Jei Je' 


i I' r 


defined in appendix 


are respectively 


\v, 


Jei Je 



(III.107) 


3 

2 


% /6(_)>-v-fe e /+1 


i W je W je , 

m 


Kl(v) 






and 


, I v;j e ,je',...-,k e ',...) = 


K 2 {v) 


^ jhii' {.CijkJv,e,jJv,e',k) j'k' Jy ,e ,j' ,e',k ') 

$klJv,e,kJv,e',l 


7r — arccos 


$kk' Jy,e,kJv,e,k' \I Jy,e' ,1' _ 


(III.108) 


1^5 jei je' i • • • i foe* i • • • ) 


^^-Vctieije'ike') 


2n 


OLv 


@{jei je' 1 k e ' ) j |u, jei je' i ■ ■ ■ j k e ' ,...), 


where 


c(jeije'i he' ) • j e (je T l)je' (Je' T 1) 


k e \k e ' + 1) - je {je + 1) ~ je'{je' + 1) 


k e '(k e ' + 1) - je (je + 1) - je' {je' + 1) 


(III.109) 


d(jeije'ike') '■= TT - arCCOS 


ke'(k e ' + 1) - jejje + 1) - je'jje' + 1) 
2\/ je{je + l)je' {je' + 1) 


(III.110) 
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Notice that while the operator H^ e e , in the Euclidean part is a graph changing operator, 

hence not preserving the original intertwiner space, the operator , in the Lorentzian part 

i s diagon al on the basis adapted to the pair of edges {e,e'}. Also, from equations (III.107) and 
(III. 108), we can deduce that the domain of the Hamiltonian operator H admits an orthogonal 
sum decomposition in terms of stable subspaces under repeated action of H. This result gener¬ 
alizes to other symmetrizations than the one proposed in (III. 105), and it may be of considerable 
importance in the elaboration of self-adjointness proofs and the calculation of the evolution of 
physical states in this model. 


IV. SUMMARY AND OUTLOOKS 

We considered a model of Einstein gravity coupled to a free scalar held, in which the dy¬ 
namics of the gravitational held is described by deparametrization with respect to the scalar 
held. In the corresponding quantum theory, constructed using the techniques of loop quantum 
gravity, the quantum dynamics is given by the evolution of the physical (i.e. gauge and dif- 
feomorphism invariant) states of the gravitational held with respect to the scalar held. This 
evolution is governed by a physical Hamiltonian operator, which we constructed in this paper. 
The implementation of the Lorentzian part of our Hamiltonian is based on regularization used 
to define the curvature operator introduced in [33]. As to the Euclidean part, we rehned and 
made precise the idea, hrst considered in m, of regularizing the curvature by means of loops 
attached to pairs of edges at a vertex of a spin network graph. 

By carefully specifying the properties of the special loops created by the Euclidean Hamilto¬ 
nian operator, we were able to define an operator which is diffeomorphism invariant, and whose 
adjoint operator is densely dehned. The second property is crucial in that it allows to sym¬ 
metrize the operator and eventually to construct self-adjoint extensions. Our regularization of 
the Euclidean term can also be applied in vacuum loop quantum gravity to define a Hamilto¬ 
nian constraint operator for which the adjoint operator is densely dehned. This question will 
be treated in a future work. On a practical level, an important feature of our Hamiltonian is 
that the volume operator does not appear in it. This implies a considerable simplification of the 
calculation of the action of the Hamiltonian on spin network states. 

The construction presented in this paper gives us a concrete and tractable Hamiltonian 
operator for loop quantum gravity coupled to a free scalar held. This makes it possible to test 
the dynamics of the theory, as the time evolution of spin networks under this Hamiltonian can 
be computed. In particular, a question of interest will be to study the evolution of semiclassical 
states describing e.g. cosmological spacetimes. 
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Appendix A: Graphical calculus for SU( 2) 

In this appendix we give the technical tools of SU (2) representation theory that are needed to 
evaluate the action of the Hamiltonian on a spin network state. For a more detailed presentation, 
we refer the reader to [39] . 
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1. The epsilon tensor 

The epsilon tensor is the fundamental invariant tensor of SU(2 ). For the spin-j representa¬ 
tion, it is given in the standard basis 6 {\j,m)} by 

e w=(-irx,-„ (A.i) 


and is represented graphically as 

(j) _ J 

e mn ~ m - > n 

It satisfies the symmetry relation e£m = (-l) 2 ^™, i.e. 

_ l _ = (-1) 2J _i 


(A.2) 


(A.3) 


The tensor gO) mn is defined to be numerically equal to emn- The contraction of two epsilons 
gives 

44'“”" = C (a.4) 

Graphically, is represented by a line with no arrow, and the above relation reads 

- > < - = - (A.5) 

Indices of SU(2) tensors can be raised and lowered using the epsilon tensor. Our convention is 
the following: 


Vm. — ^ 


(A.6) 


2. Intertwiners 

The intertwiner between three representations j \, j 2 and js is given by the Wigner 3j-symbol: 



(A.7) 


() All the relations and conventions given in this appendix refer to the basis {\j, rn)}, which is the usual eigenbasis 
of the operators J 2 and J 3 . For j = the basis {| |, +|) , 11, — |)} = (|+) , |—}} of the space Hi /2 is defined 
by specifying that T 3 G su( 2 ) and e £ ® H*/ 2 are represented by the following matrices: 


t \A _ * 

(73) B ~ ~ ^ 


1 0 
0 -1 


cab = 


0 1 
-1 0 


the indices A and B taking the values + and —. For a general representation j, the normalized states | j,m), 
which span the space Hj, are obtained as symmetric tensor products of the states |+) and |—): 


| j, m) = J ^ + ^ 2 ^) 1 —~“ (l+) ® ® 1 +) ® |-) ® ® |-} + all permutations^. 

j + m states j — m states 
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The order of the spins in the symbol is indicated by a + or a — at the node. Thus, 



(A.8) 


which is a graphical representation of the relation that switching two columns in the symbol 
(A.7) multiplies the symbol by (— \y i+ te+]3_ Another symmetry relation is 




(A.9) 


When one of the spins is zero, the 3j-symbol reduces to the epsilon tensor: 



(A-10) 


where dj = 2j + 1 is the dimension of the representation j . Intertwiners of higher valence are 
constructed by contracting several three-valent intertwiners. For example, the objects 


,(k) 

mi 7712 5 m 3 7714 


f h 32 73 J4 \ 

4 -^ \rri 1 m2 /r ) /j, m3 m^J 



(A.ll) 


form a basis in the space of intertwiners between the representations j\. j'2, j'3 and j’4. (Note 
that Lm\m 2 -,m 3 m 4 is not normalized; its norm is 1 /\/~dk-) 

In calculating the action of the Hamiltonian, a specific basis of intertwiner states at the vertex 
v is chosen. These states are denoted by | v; j e , j e ', je" , ■ ■ - k e ", ...) and defined as 


jei je' j je" i ■ • ■ i k e' > k e ", . . . ) 


je 



+ 


+ 


(A.12) 
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3. 6 j- and 9j-symbols 


A contraction of four 3j-symbols defines the 67 -symbol. 



(A.13) 


This object appears 
and which gives the 


in the following relation, which is a special case of equation (A. 19) below, 


change of basis between two different bases of the form (A. 11 ): 




£>(-!)*+» 


— k—x ] 7I 72 

74 73 


72 73 

Similarly, the 9j-symbol is defined as 




7i 72 
ki k 2 

h h 




It can be expressed in terms of 67 -symbols by the relation 


7i 72 
k\ k 2 

h h 





(j 1 72 73 1 

\h h X j 



k 2 

x 


h\jh 

h J \x 


h 13 1 

7i h J ' 


(A.14) 


(A.15) 


(A.16) 


4. Expanding invariant tensors 

Any invariant tensor 7 mi ... mjv , having indices in representations j\ . ■ ■ ■ ,jN, is an element of 
the space of intertwiners between the representations j 1 ,..., j'at, and as such, it can be expanded 
using a basis of the intertwiner space. Expressing a tensor with N indices as a block to which 
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N lines are attached, one has the relations 



(A-17) 


(A.18) 


(A-19) 


as well as the straightforward generalization of the last relation for tensors of higher order. 


5. Group elements 


The representation matrix for a group element is expressed graphically as 


D^ m n (g) 


n 




m 


(A.20) 


The inverse matrix is given by the relation 

D^nig- 1 ) = e^e^D^(g), (A.21) 


or 








(A.22) 


In computing the action of the Hamiltonian, we need to know the action of a flux operator 
Ps,i on a holonomy h e . In the cases where the intersection v between the surface S and the edge 
e is the beginning or ending point of the edge, this action is given by 


Ps,iD^ m n (h e ) = S,e)J v , e ,iD^ m n (h e ) 


= ^(S,e) x 




if e begins from v 
if e ends on v 


(A.23) 


where n(S, e) = +1 if the orientation of e agrees with the direction of the normal vector of S, 
and n(S,e) = — 1 if the orientations of e and S are opposite. 

( j) 

The matrix r- is proportional to an intertwiner between the representations j, j and 1; the 
precise relation is 
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i ' n 


Wj n 


J + J 

- • > 


where Wj = iy/j(j + 1 )(2j + 1). Therefore we can write (A.23) graphically as 

j : J 

2 


J 


Ps, l -h>- = ^W jK (S,e) V — -• > 


1 


Ps,i 


D> 




D> 


+ 


v 

-* -• 


(A.24) 


(A.25) 
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